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(54) Title: METHOD FOR MODELLING A SURFACE AND DEVICE FOR IMPLEMENTING SAME 
(54)Titre: PROCEDE DE MODELISATION D'UNE SURFACE ET DISPOSITIF POUR SA MISE EN (EUVRE 
(57)Abstract 

A method for obtaining a model of a surface, characterized by the steps which involve obtaining measurements of geomet- 
rical data concerning specific points on the surface, making a grid of the surface, with the grid passing through said points, mem- 
orizing, at an address which is specific to each node of the grid, the coordinates of the node, the number of satellites of the node, 
information for access to the addresses of said satellites and thereafter to information which relates to them, and geometrical data 
which may be associated with said node; for each node, defming a local roughness index obtained from a weighted sum of the 
current coordinates of the node and its satellites, defining the sum of an overall roughness index representing the sum of all the 
local roughness indices, and of an overall mdex of the infringement of said geometrical data, iteratively adjusting the coordinates 
of indefinite nodes, by using at each adjustment the sum of a weighted combination of current node neighbour coordinates and of 
a combination of geometrical data associated with the node, in order to minimize said sum, and creating a mode! of the surface 
on the basis of the adjusted coordinates. 

(57) Abrege 

L'invention conceme un procede d'obtention d'une modelisation d'une surface. II est caracterise par les etapes consistant a 
obtenir des mesures de donnees geometriques relatives a des points specifiques de la surface, efTeauer un maiilage de la surface 
passant par lesdits points, memoriser a une adresse specifique a chaque nceud du maiilage les coordonnecs du naud, le nombre 
de satellites du nceud, une information d'acces aux adresses desdits satellites et par suite aux informations y relatives et des don- 
nees geometriques eventuellement associees audit nceud; pour chaque naud, deflnir un indice de rugosite local obtenu a partir 
d'une somme ponderee des coordonnees actuelles du ncEud et de ses satellites, deflnir la somme d'un indice de rugosite global, 
constitue par ['addition de tous les indices de rugosite locaux, et d'un indice global de violation desdites donnees geometriques, 
ajuster de fapon iterative les coordonnees des nceuds indetermines, en utilisant a chaque fois I'addition d'une combinaison ponde- 
ree des coordonnees actuelles de voisins du nceud et d'une combinaison des donnees geometriques associees au nceud, pour ame- 
ner ladite somme a un minimum, et creer une representation de la surface a partir des coordonnees ajustees. 
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PROCEDE DE MODELISATION D'UNE SURFACE ET DIS POSITIF POUR 
SA MISE EN OEUVRE . 

La prfesente invention concerne d'une fa<?on 
gen^rale 1 • investigation dans des corps k trois 
dimensions, et concerne en particulier un nouveau proc6d6 
pour I'obtention d'une representation d'une surface 
situ6e au sein de ce corpse partir d'un ensemble limits 
de donn^es g6oni§triques connues relatives h ladite 
surface. 

L' investigation dans des corps tridimensionnels 
est une preoccupation majeure tout particuli^rement dans 
les domaines de la g6ologie, de la g6ophysique et de la 
biologie. 

Ainsi, par exemple dans le domaine de la 
g6ophvsigue, on souhaite obtenir des representations 
aussi fiddles que possibles de surfaces situ6es par 
exemple h 1' interface entre deux zones de natures ou de 
propri6t6s diff6rentes, k partir de donn6es obtenues lors 
de campagnes de prospection. voire en cours 
d' exploitation d'une richesse souterraine. 

En effet, l'efficacit6 de 1 • investigation du corps 

20 tridimensionnel. et en particulier de la recherche de 
p6trole ou autre richesse souterraine. est conditionn^e 
par la fid61it§ avec laquelle on peut reconstituer et 
repr6senter ce genre de surface. 

Par ailleurs, dans , le domaine medical et 

25 biologique, I'on sait actuellement obtenir, par des 
proc§d6s divers, des representations de coupes d'un corps 
vivant ou analogue, et I'on recherche 6galement 
reconstituer et ^ repr^senter, de facon 
tridimensionnelle, une surface situ^e ^ 1' interface entre 

30 deux milieux, et notamment les contours d'un organe ou 

autre. , 

II existe k I'heure actuelle un certain nombre de 
techniques de mod61isation permettant d' obtenir la 
representation d'une surface dans un corps i troxs 
35 dimensions que I'on cherche k exploiter ou^ traiter. On 
citera en particulier les interpolations de Be^ier.ou les 
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fonctions splines {on se r6f^rera pour plus de details 
par example ^ 1 ' ouvrage Geometric Modelling, M, E- 
Mortenson, Ed. John Wiley, 1985) . 

Cependant, ces techniques connues sont toutes tr^s 
5 mal adapt^es ^ la prise en compte de donn^es h6t6rog6nes , 
^ savoir des donn6es relatives k la fois aux coordonn^es 
de points de la surface, k des orientations imposfees de 
la surface, ^ des liens g6om6triques entre deux points 
distants, etc. En outre, ces techniques sont inadaptfees k 

10 

une trop grande dispersion des surfaces ou encore & des 
anomalies de surfaces telles que des plis , des cassures , 
ou encore des surfaces discontinues. Plus pr6cis6ment, 
dans certains cas , les fonctions mises en oeuvre dans ces 
techniques soit ne convergent pas, soit n'admettent 

15 aucune solution, soit admettent au contraire une 
plurality de solutions. 

Enfin les techniques connues n' off rent pas la 
possibility de prendre en compte la notion de plus ou 
moins grande certitude avec laquelle certaines donn^es de 

20 la surface sont connues* 

La pr6sente invention vise k pallier les 
inconv6nients des proc6d6s de la technique ant6rieure et 
^ proposer un proc6dfe qui puisse prendre en compte des 
informations g6om6triques non homog^nes et/ou trfes 

25 dispers6es, ainsi que les anomalies pr6cit6es et 
d'autres, en off rant ^ chaque fois une mod61isation 
unique et sans Equivoque de la surface. 

Un autre objet de 1" invention est de proposer un 
proc6d6 qui puisse s ' appuyer sur des donn6es g^om^triques 

30 prenant en compte simultan6ment un degr§ de certitude 

ou de precision li6 ^ ces donn6es. 

A cet effet, la pr6sente invention concerne tout 
d'abord un proc6d6 d'obtention d'une mod^lisation d'une 
surface representative notamment de 1' interface entre 

2j deux zones de natures ou de propri6t6s diffferentes dans 
un corps k trois dimensions tel qu'une. formation 
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g6ologique ou un corps vivant, caract6ris6 en ce qu'il 
comprend les 6tapes consistant k: ' 

- obtenir au moyen d'appareils de mesure un 
ensemble de donn6es g6oin6triques relatives k la surface 

5 et associ6es respectivement k des points donn6s de ladite 
surface; 

- effectuer un maillage de la surface de telle 
sorte que 1' ensemble desdits points donn6s soit un sous- 
ensemble des noeuds du maillage; 

10 _ m^moriser h une adresse donn6e d'une ni6inoire , 

specif ique h chaque noeud du maillage, les informations 

suivantes : 

* les coordonn6es du noeud consid6r6, 

* le nombre de noeuds satellites du noeud 

15 consider 6, 

* une information permettant d ' acc6der aux 

adresses sp6cifiques desdits noeuds satellites et par 
suite aux informations relatives h ces derniers , 

* le cas 6ch6ant des donn6es g6om6triques 
associfees audit noeud consid6r6, 

- pour chaque noeud du maillage, d^finir un indice 
de rugosity local obtenu h partir d'une somme pond§r6e 
des coordonn6es actuelles du noeud et de ses satellites, 

- d^finir la somme d'un indice de rugosit6 global, 
constitu6 par 1' addition des indices de rugosity locaux 
associ6s h chaque noeud, et d'un indice global de 
violation desdites donn6es g^omfetriques ; 

- ajuster de facon iterative les coordonn6es de 
chacun des noeuds dont on ne connalt pas les coordonn6es 
pr6cises, en utilisant ^ chaque fois 1' addition d'une 
combinaison pond6r6e des coordonn6es actuelles des 
satellites et des satellites des satellites dudit noeud 
et d'une combinaison des donn6es g6om6triques associ6es 
audit noeud, de mani4>re h amener ladite somme k un 

32 minimum ; 
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- cr6er une representation de la surface h partir 
des coordonn6es ajust^es de chacun des noeuds . 

L* invention concerne 6galement un dispositif pour 
I'obtention d'une modelisation d'une surface 
5 representative notamment de 1' interface entre deux zones 
de natures ou de propri6t6s diff6rentes dans un corps ^ 
trois dimensions tel qu'une formation g6ologique ou un 
corps d'etre vivant, caract6ris6 en ce qu'il comprend : 

- des moyens de mesure un ensemble de donn6es 
g6ometriques relatives A la surface et associ^es 
respectivement ^ des points donnas de ladite surface; 

- des moyens pour effectuer un maillage de la 
surface de telle sorte que 1' ensemble desdits points 
donnfes soit un sous-ensemble des noeuds du maillage; 

15 - une m^moire contenant ^ des adresses donn^es , 

sp6cifiques ^ chaque noeud du maillage, les informations 
suivantes : 

* les coordonnees du noeud consider^, 

* le nombre de noeuds satellites du noeud 
20 considfere, 

* une information permettant d'acc6der aux 
adresses sp6cifiques desdits noeuds satellites et par 
suite aux informations relatives k ces derniers, 

* le cas ech6ant des donn6es g^ometriques 
25 associ6es audit noeud consid6r6, 

- des moyens de calcul pour ajuster de fagon 
iterative les coordonnees de chacun des noeuds dont on ne 
connalt pas les coordonnees pr6cises, en utilisant k 
chaque fois 1* addition d'une combinaison ponderee des 

2Q coordonnees actuelles des satellites et des satellites 
des satellites dudit noeud et d'une combinaison des 
donnees geometriques associees audit noeud, de maniere k 
amener k un minimum une somme d'un indice de rugosite 
global f constitue par 1 ' addition d ' indices de rugosite 

22 locaux respectivement associes i chacun des noeuds et 
chacun obtenu k partir d'une somme pond6ree des 
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coordonn6es actuelles dudit noeud et de ses satellites, 
et d'un indice global de violation desdites donn6es 

g6oin6trigues ; et 

- des moyens pour cr6er une representation de la 
5 surface k partir des coordonn6es ajust6es de chacun des 
noeuds . 

D'autres aspects,' buts et avantages de la pr6sente 
invention apparaitront mieux la lecture de la 

description dfetaillfee suivante d'une forme de realisation 
10 pr6f6rfee de celle-ci, donn^e k titre d'exemple non 
limitatif et faite en r6f6rence aux dessins annexes, sur 
lesquels : 

la figure 1 est un exemple de graphe incorporant 
un maillage utilise pour modfeliser une surface, ce graphe 
etant defini par 1' ensemble de tous les c6t6s de facettes 
dont les sommets constituent 1* ensemble Q d^crit plus 
loin, et illustre 6galement des f Inches symbolisant des 
deplacements de noeuds du maillage effectufes 
interactivement par 1 ' utilisateur pour faconner la 
surface; 

_les_figure 2(a) & 2(d) illustrent quatre cas de 
contr61e de la forme d'une surface t I'aide de quatre 
types differents de contraintes vectorielles; les 
vecteurs A^u sent supposes connus directement sur les 
figures 2(a) k 2(c) et connus par leur orthogonality h un 
vecteur V donn6 sur la figure 2(d); 

la figure 3 illustre un exemple de surface 
geologique complexe (d6me de sel) inod61is6e k I'aide du 
precede de la pr6sente invention; 

la figure 4 illustre un exemple de surface 
^° biologique complexe (cerveau d'embryon) mod61isee h 
I'aide du procede de la presente invention k partir de 
donnees ponctuelles obtenues par une succession de 
coupes ; 

la figure 5 illustre 1 ' organisation de la m6moxre 
associee k un atome A(k) de noyau ; 



20 



25 



35 
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la figure 6 illustre un premier exemple de 
memorisation d'une orbite d'un atome donn6, dans un 
tableau de nic cases; 

la figure 7 illustre un second exemple de 
memorisation d'une orbite d'un atome donne, dans nk 
tableaux chaln6s de deux cases chacun; 

la figure 8 illustre la memorisation d'une donn^e 
g6om6trique ou contrainte; 

la figure 9 illustre la mani^re de m6moriser un 
ensemble de donn6es g6ometriques ou contraintes relatives 
h un atome A(k) donn^; 

la figure 10 illustre la memorisation dans un 
tableau de 1' ensemble des atomes relatifs A une surface 
S; et 

1^ figure 11 illustre 1 ' enchalnement de 1 ' acc^s 
audit ensemble d' atomes. 

On indiquera tout d'abord que, d'une figure ^ 
1" autre, des elements ou parties identiques ou similaires 
sent designes par les mSmes num6ros de reference. 
20 En outre, on va indiquer ci-dessous un certain 

nombre de publications dont les contenus respectifs sont 
k considerer comme inclus dans la pr6sente description. 

(1) A method of Bivariate interpolation and smooth 
surface fitting based on local procedures, H. AKIMA, J. 

25 ACM 17.1, 1974; 

(2) Shape reconstruction from planar cross-sec- 
tions, J.D. BOISSONNAT, ACM Transactions on Graphics, 
vol. No. 2, 1986; 

C3) Machine contouring using minimum curvature, 
30 I-C. BRIGGS, Geophysics 17,1, 1974; 

(4) Triangular Bernstein-Bezier patches, G. FARIN, 
Computer-aided geometric design, 1986; 

(5) Primitives for the manipulation of general 
subdivisions and the computation of Voronoi diagrams, L. 

35 GUIBAS et J. STOLFI , ACM Transactions on Graphis , vol. 
No. 2, 1985; 
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(6) Three-dimens-ional graphic display of disconnected 
bodies, J.L. MALLET, Mathematical Geology, vol. 20, no. 8, 
1988; 

(7) Geometrical modelling and geostatistics , J.L. 
5 MALLET, Proceedings of the Third International 
Geostatistics Congress,' Kluwer Academic Publishers, 1989; 

(8) Discrete Smooth Interpolation, J.L. MALLET, ACM 
Transactions on Graphics, April issue, 1989; 

(9) Geometric Modelling. M.E. MORTEN SON . John Wiley, 

10 1989. 

Par ailleurs, la description qui suit d ' une forme de 
realisation pr^fer^e de 1 • invention est elle-meme suivie 
d'une annexe exposant une m6thode de reconstitution d'une 
surface d6riv6e de la m6thode expos6e dans le document (8) 
15 et sur laguelle est bas6 le proofed^ de la pr6sente 
invention. 

1. DESCRIPTION DU PROCEDE DE CODAGE ET D'AJUSTEMENT 

DE SURFACES COMPLEXES SELON L' INVENTION 
1.1. Introduction 
20 1.1.1 Ob jet de la description 

on va d6crire ci-dessous un proc6d6 industriel de 
codage et d'ajustement de surfaces complexes rencontrfees 
par exemple lors de la inod61isation d' interfaces entre 
zones de natures ou de propri6t6s diff^rentes de couches 
25 g^ologiques ou de corps d'organismes biologiques. 

Le r^sultat de ce codage peut §tre transcrit sous forme 
d- images et/ou de maquettes en matferiel solide (par exemple 
plastique) k 1 ' aide d ' un ordinateur graphique ou de tout 
autre proc6d6. En g^ologie. par exemple, ces images et/ou 
30 maquettes sont ensuite utilis6es pour optimiser 
1- exploration, 1 'exploitation et la gestion des ressources 
du sous -sol telles que : 
. pfetrole 
. mineral 

35 • 

. etc. 
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Le proc^d^ d^crit peut fitre decompose en deux sous- 
procedes complementaires : 

1. un sous-procede de codage d'une surface 

2. un sous-procede d'ajustement d'une surface aux donnees disponibles 
hn geologie, par exemple, ces donnees peuvent etre : 

• la position (exacte ou approchee) de points sur une coupe de la 
surface par un plan 

• la position (exacte on approchee) de I'intersection d'un forage avec 
la surface geologiquecorrespondant a la separation entre deux couches 

• le plan tangent (exact ou approche) en un point donne de la surface 
correspondant a la sq>aration entre deux couches 

• le rejet (exact ou approche) d'une faille (=cassure) affectant la 
surface correspondant a la separation entre deux couches 

• etc ... 

du nrocedTt'"^' codage est con^u de fagon a optixniser les performances 
fon^ ' f ajustement base sur une variante de la" methode DSI (cf ref 
[8]) specialement adaptee a la inodeUsation geometrique des surfaces. Cette 
vanante de la methode DSI utiKsee ici n'a encore jamais ete pubUee et ^ t 
presentee en annexe. 



10 



20 1.1.2 Decomposition d'une surface S en facettes 

Soit 5 une surface a modeliser. Comme le montre la Figure 4 et ainsi qu'il 
a ete suggere dans les references [6] et [7], celle-ci sera supposee approchee 
par un ensemble de facettes triangulaires et/ou polygonales. Ces facettes sont 
caractensees pax : 

25 •La osnnaissance des coordonnees de leurs sommets constituant un en- 
semble fini de N points „ . . . , <pj^} reguUerement repartis sur S. Ces 
sonmiets sont numerotes de 1 jusqu'a N et les coordonnees du fc^me 
sommet seront notees (Vfc,Vfc.Vfc) 

• La donnee de leurs cotes constitues par des Uens entre couples de som- 
30 mets {<^i,<pj) 

L'ensemble des cotes des triangles et/ou polygenes composant une surface 
S constitue un graphe G dont les noeuds sont precisement les sommets des 
triangles et/ou polygenes. Pour cette raison, et afin d'eviter des periphrases 

FEUILLE DE REMPLACEMENT 
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2) <?,-€A(fc) 



inu tiles, nous utiliserons de fa?on equivalente le mot "noeud" a la place de 
"sommet de triangle et/ou polygene" : 

noeud = { sommet de triangle et/ou polygene } 

1.1.3 Definitions preliminaires 

5 Ensemble des satellites A(/:) 

Soit (pi, le k^^^ noeud constituant une surface 5. Nous appellerons ensemble 
des satellites de (pt, i^ous noterons A(fc) I'ensemble des noeuds differents 
de tels que ^t) soit un cote de triangle ou de polygene. En d'autres 
termes, I'ensemble A(A:) est tel que : 

1) A(fc) = {<^^,,...,<?,„J 

(pj et extremites d'un 

< meme cote de triangle ou de poly- 
. gone 

Comme indique dans la relation ci-dessus, nous conviendrons de noter Uk le 
15 aombre de satellites r attaches au noeud ip^. 

— Cette notion d'orbite introduite dans les articles [6] et [7] n'a jusqu'a 
present pas fait Tobjet de precede de codage particulier; une partie imper- 
tante de 1 ' invention consists h proposer un proc6d6 de 

codage pour A(A:). 

70 Voisinage iV(fc) 

Soit le k^^^ noeud d'une surface 5. Nous appellerons voisinage de (p^ et 
nous noterons N{k) I'ensemble des noeuds constitues par la reunion de (p^ et 
de ses satellites. En d'autres termes, on a : 

^(fc) = U A(fc) 

25 Nous appellerons voisin du noeud tout noeud appartenant a I'ensemble 
Ar(fc). En d'autres termes, on a : 

€ N{k) <^ {(pj = Voisin de ip^} 



Cette notion de voisinage introduite dans les articles [7] et [8] n'a jusqu'a 
present pas fait I'objet de precede de codage particulier; une partie importante 
30 de 1 • invention consiste k proposer un proc6d6 de <?(>dage 
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Atome A{k) et noyau associe 

Soit le **me ^^^^ ^.^^^ ^^^^ ^ ^^^^ appellerons atome de noyau 
♦ du noeud <pf. 

t derensembleA*(A:)desadressesotidescodes|/l*{'i,i ' A*(-i W r.^. 
metta^t de retrouver da^ une memoire I'ensemble d^s atorkes {A(l C 
A (Jnjicorrespondants dont les noyaux sont les satellites de 
En d'autres termes, on a : 

A{k) = (^.,A*(A:)) = (^,,{A*(iO,...,A*OnJ}) 

rik = nombre de satellites de <^^, 

A*{j) = Adresse ou code d'acces de Tatome 
A{j) correspondant au noyau du j^"^® 
satellite de 

Pot!^tT'°'' "^if^T I'^ticle [6] n'a jusqu'a present pa^ fait 

objet de procede de codage particuUer; une paxtie importante de 
1 invention consiste ^ proposer un proc6d6 de codage pour A(fc) 

Contraintes et types de contraintes 

Nous appeilerons contrainte relative a un atome A{k), toute condition relative 
a la position du noeud correspondant au noyau de cet atome. Par exemple 
en geologic on pent citer les contraintes suivantes qui sont d'une importance 
capitale pour I'ajustement dW surface aux donnees sures ou approchees 
observees ou mesurees par des procedes manuels ou geophysiques : 

• Les ooordonnees {^^1,^1,^%) du noeud (^^^ sont connues de fagon exacte 
Une teUe contrainte est dite de type noeud de controls et est utilisee 
en geologic par exemple pour coder I'intersection de la surface modelisee 
avec un forage. 

• Les coordonnees du noeud sont connues de facon ap- 
prochee avec un coefficient de certitude donne. Une telle contrainte 
dite de type noeud de controle flou et est utiUsee pour coder la posi- 
tion approximative d'un noeud de la surface modelisee supposee mesuree 
de fagon approximative par un procede donne. Par exemple, 
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— en geologie, les donnees sismiques sont de ce type, 

— en biologie, les donnees de coupes d'organes observees sous micro- 
scope ou par echographie ou encore par scanner sont de ce type. 

• Le vecteur — tp^) joignant le noeud <^jt ^ autre noeud est 
connu de fagon approchee avec un coefficient de certitude donne. Une 
telle contrainte est dite de type lien vectoriel flou et est utilisee en 
geologie par exemple pour coder le rejet d'une faille, le noeud ^^etant 
sur un bord de la faille et le noeud (px ^tant alors sur le bord oppose 
(voir Figure 2). 

• Le vecteur V orthogonal au plan tangent a la surface S au point corre- 
spondant au noeud c^j^ connu de fagon approchee avec un coefficient 
de certitude donne. Une telle contrainte est dite de type normale floue 
et est utilisee en geologie pax exemple pour coder le pendage des couches 
mesure par des procedes manuels ou geophysiques (pendagemetrie). 

15 # etc ... 

Les coeflBcients de certitude mentionnes ci-dessus sont des nombres positifs 
dont la valeur est censee representer la confiance que Ton peut avoir "a priori'' 
dans les donnees qui lui sont attachees. La prise en compte des contraintes 
floues pr6sent6es ci-dessus constitue l*une des 
20 caract§ristiques de 1' invention est pr^sent^e en detail 
dans 1 ' annexe (voir annexe - section 2.5 et 2.7.6.). 

Con vent ron 

Afin de simplifier les notations, on convient d 'identifier les numeros des noeuds 
d'une surface 5 avec les noeuds eux-memes. II s'ensuit que Ton a : 

25 # = ensemble des satellites -{(^j^, ,<^j-^^} de ipf^ ou ensemble des 

numeros {ji, . . . , jnt} correspondants 

• N{k) = ensemble des voisins {^j^^^^ii? - - ^jnj^} ou ensemble des 
numeros {A:, ji, . . . , jnk} correspondants 



5 



10 



Etat de Tart 

30 Les methodes classiques 

On trouvera dans la reference [7] une description tres complete de Tetat de Tart 
concernant la modelisation de surfaces complexes. En particulier, dans cet— 
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article, on expose les raisons pour lesqueUes les methodes claasiques utilisees 
en cartographie automatique et en Conception Assistee par Ordinateur (CAO) 
ne sont pas adaptees a la modelisation de telles surfaces comme on aurait 
pu le croire a priori. A la fin de cet article [7], on introduit brievement les 
fondements de la methode DSI (quelquefois appelee DSI/DS A) permettant de 
surmonter les difficultes inherantes aux methodes classiques. 



La methode DSI 



Cette methode est basee sur la notion de "rugosite" d'une surface S composee 
de xacettes triangulaires et/ou polygonales. Tout d'abord, la notion de rugosite 
1 0 R(jp\k) au noeud d'une surface S est definie a I'aide de la formule suivante 
(Voir ref. [7], [8] et Annexe) : 



Rim =11 E 



Les coefficient positifs ou nuls {t;-(Ar)} intervenant dans la formule ci-dessus 
sont d^ ponderateurs choisis par I'utilisateur charge de modeliser la sur- 
face Paxmi I'infinite de choix possibles, I'un des plus interessants appelle 
ponderation harmonique" consiste a poser : 



si a e A(fc) 

Les rugosites locales ainsi definies sont ensuite combinees en une rugosite 
20 globale R{<^ definie conmie suit : 

k 

Les coefficients positifs ou nuls fM{k) intervenant dans la formule ci-dessus 
servent a moduler la contribution des rugosites locales i2((^jfc) a la rugosite 
globale R{ip). Si i'on decide par exemple d'utiliser une ponderation uniforme 
25 il suffit de poser : ' 

Parmi les autres choix possibles pour les coefficients ^(jb), on peut decider 
d'utiUser pax exemple les ponderateurs suivants ou m est une constante donnee : 



fi{k) = \ ^ '?k est un point de donnee 

\ m > 1 si <^;t II 'est pas un point de donnee 
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Le principe de la m^thode DSI consists ^ ajuster la position 
de chacun des noeuds 1^ surface S* iaod^lis6e de telle 

sorte que celle-ci soit la plus lisse possible tout en res- 
pectant les donn^es et contraintes r^gissant sa forme. 
5 article [8] mentionn^ en bibliographie est consacr^ ^ une 

presentation math^matique approfondie de la m^thode DSI et 
propose une technique de minimisation de la rugosity globale 
respectant les donn^es et contraintes. 

Version interactive de la nnethode DSI 

10 La methode de minimisation de la rugosite globale de la methode DSI proposee 
dans rarticle [8] est mal adaptee a une utilisation interactive sur calculateur 
electronique (ordinateur ou microprocesseur specialise), c'est pourquoi nous 
proposons en annexe une autre technique de nunimisation originale sur laque- 
lle est basee 1' invention . Ainsi qu'il a ete dit precedenmient, Tobjet de 

1 5 1" invention est precisement de proposer un procede de codage et un procede 
d'ajustement de la surface 5 suivant cette nouvelle version de la methode DSL 

1.2 Exemples d'images d'objets 

codes et ajustes par notre procede 

1.2.1 Exemple d'image de surface geologique 
2 0 olStenue par notre procede 

La recherche du petrole est generalement effectuee de fa^on indirecte en ce sens 
que plutot que de rechercher le petrole, on cherche des couches geologiques 
reservoir suceptibles de le contenir ou de le capturer; par exemple, la deformation 
des couches de sel donne naissance a ce que Ton appelle des "domes de sel" qui 

2 5 sont consideres conmie generant d'excellents pieges a petrole. Pour pouvoir 
localiser avec precision la position du petrole il est necessaire de connaitre 
la forme exacte des couches "pieges'' et malheureusement dans de nombreux 
cas il n'existait jusqu'a ce jour aucun procede pratique permettant de les 
modeliser de fagon satisfaisante; le cas des domes de sel est a ce titre considere 

30 comme Tun des plus complexes et pour prouver I'efEcacite de notre procede 
de codage et ajustement nous presentons sur la Figure 3 Timage d'une telle 
surface obtenue avec celui-ci. 
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1.2.2 Exemple d'image de surface biologique 
obtenue par notre procede 

Les procedes d'acquisition de donnees tels que scanner ou echographie utilises 
en medecme sont mathematiquement ident4ques aux procedes d'acquisition de 
donnees sismiques utilises en prospection petroliere. Par exemple, Techographie 
est identique a la sismique reflexion et le fonctionnement des scanners est 
base sur la tomographie tout comme la tomographic sismique Des lors il 
n est pas etonnant que notre procede de codage et ajustement des surfaces 
soit egalement utilisable pour modeliser la peau d'organes biologiques et nous 
proposons sur la Figure 4 le modele du cerveau d'un embryon humain de 7 mil- 
limetres de long obtenu a I'aide de notre procede de codage et d'ajustement 
II est h noter qu'un aussi petit organe (le cerveau fait 
1 mxlli„.^tre de long) ne pouvait jusqu'^ present gtre 
observe que sous forme de coupes microscopiques 



20 



1.3 Procede de codage d'lme surface S 

1.3.1 Notion de memoire et d'adresse de memoire 

Nous appellerons "memoire" tout dispositif de calculatrice electronique qui 
enregistre des donnees, des codes ou resultats partiels qu'il restitue au mo- 
ment opportun pour les faire intervenir dans la suite des operations De plus 
nous appellerons "adresse" d'une memoire, tout code permettant de retrouver 
1 emplacement physique de la memoire a I'interieur de I'ordinateur. 

1.3.2 Codage d'un atome A(k) 

Nous proposons de coder et stocker chaque atome A(ib) dans un ensemble de 
25 memoires a Taide du procede suivant ou I'une de ses variantes : 

Procede de codage 

Conome le suggere la Figure 5, dans une zone memoire contigue commengant 
a radresse A*{k), on stocke de fa^on consecutive les informations suivantes 
relatives a I'atome A{k) correspondant au noeud : 

30 1- noyau = 
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• soit les coordonnees {<pty^t^^k) noeud (pf^ 

• soit I'adresse ou le code ou toute autre information permettant de 
retrouver les coordonnees {^l^^l^^k) noeud (pf^ 

2. nb«sat = 

5 • soit le nombre rik de satellites lies du noeud 

• soit I'adresse ou le code ou toute autre information permettant de 
retrouver le nombre de satellites lies du noeud (pj^ 

3. sat = 

Adresse d'une zone memoire contenant les adresses ou les codes d'acces 
1 0 {^*0i )^ • • - ^ -^*Onfc )} permettant de retrouver les rik atomes {A{ji ),..., 

^ OnJ}^^^^ noyaux {^j,, - - y^j^^} sont les satellites du noeud tp,^ 

4. const = 

Codage de I'ensemble des contraintes (cf section 1.1,3) pouvant affecter 
le noeud t^j^. Ce codage sera presente en details dans la section 1.3.3. 

15 5. info = 

Zone memoire de taille variable contenant zero, une ou plusieurs infor- 
mations complementaires concemant I'atome A{k). Ces informations 
sont specifiques a I'application particuliere envisagee; par exemple en 
prospection petroliere, si A(^k) est un atome correspondant a la surface 
2 0 separant deux couches Ci et Cj, alors il peut etre interessant de stocker 
dans le champ info des informations telles que : 

• liste des proprietes physiques (porosites, vitesse des on des sismiques, 
...) de la couche Ct au noeud (p^ correspondant a I'atome A{k) 

• liste des proprietes physiques (porosites, vitesse des ondes sismiques, 
2 5 .-.) de la couche C2 au noeud <p^. correspondant a Tatome A(k) 

• liste des proprietes geologiques (facies geologique^ presence de petroie, 
..,) de la couche Ci au noeud (pf^ correspondant a Tatome A{k) 

• liste des proprietes geologiques (facies geologique, presence de petroie, 
...) de la couche C2 a.u noeud tpf^ correspondant a Tatome A{k) 

30 « etc - 
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Variante 1 

La zone memoire dans iaquelle sont stockees les adresses on les codes d'acces 
\^ U I A ( )j permettant de retrouver les njt atomes 4A(u) A(i U 

5 ttr?eZ™ns^ 

1. La premiere fagon illustree par la Figure 6 consists a utiliser un tableau 
de Uk memoires consecutives contenant les adresses . . . , A*U^ )} 

2. La seconde fagon iUustree par la Figure 7 consists a utiliser n, tableaux 
comportant chacun 2 memoires consecutives telles que oour le ^eme 

1 0 ces tableaux : m p ur le q de 

• la premiere memoire contienne I'adresse A*ij„) de I'atome corre- 
spondant au or®™® satellite de I'atome A{k) 

• la seconde memoire contienne I'adresse du (o + 1 )enie tableau Si a 
est egal au nombre de satellites (iib_sat=n,), alors cette seconde 
memoire est soit inutilisee, soit contient un code indiquant que I'on 
a attaint le demier satellite de A(k) 



15 



Variante 2 



II est possible de ne pa^ stocker le champ nb^at mais dans ce cas, si la liste des 

20 Ltr^'r-^^ ^^"^ .If (^■^)'- • • '^*OnJ} est codee sous forme d'un seul 
2 0 tableau (voir vanante 1 ), il faut rajouter un emplacement memoire A*(n, + 1 ) 
dans lequel on^lace un code indiquant que la liste A*(/b) est terminee. 



Variante 3 



Le champ sat pent etre choisi suffisamment grand de fa^on a y stocker di- 
rectement les adresses ou codes d'acces . . . , A*(j„^)} 



2 5 Variante 4 



Un grand nombre de variantes equivalentes sont possibles suivant la taiUe 
du champ info et de la fagon dont on redecoupe celui-ci. II est toutefois 
bien entendu que ces variantes n'ont aucune importance pour I'utilisation 
de 1 atome A{k) en connection avec I'algorithme d'ajustement base sur la 
3 0 methode DSI decrite en annexe. 
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Variante 5 

On pent rajouter d'autres champs a Tatome A{k) decrit ci-dessus mais ceux-ci 
n'ont aucune importance pour I'utilisation de I'atome A{k) en connection avec 
Talgorithme d'ajustement baise sur la methode DSI decrit en annexe, 

5 Variante 6 

Quelles que soient les variantes, la seule chose importante pour Timplementation 
efficace de Talgorithme DSI est la suivante : 

Ayant Tadresse ou le code d'acces A*{k) d'un atome il est 

necessaire d'acceder de la fagon la plus directe possible aux atomes 
1 0 dont les noyaux sont les satellites du noyau de A{k), 

L'object essentiel du procede de codage de A{k) dans une memoire decrit 
ci-dessus ainsi que toutes ses variantes est precisement de permettre Tacces 
direct aux satellites. 

1*3*3 Codage des contraintes 

1 5 Chaqu6 contrainte relative a un atome A{k) sera stockee dans une zone memoire 

specifique dont la description est proposee dans ce q\ii suit : 

Procede de stockage d'une contrainte 

Le codage d'une contrainte particuliere relative a un atome donne est efFectue 
dans ime zone memoire comprenant les deux champs contrainte et next 

2 0 comme I'indique la Figure 8 • Ces deux champs sont utilises de la fagon suiv- 

ante : 

• contrainte = 

Zone memoire contenant : 

— soit les informations relatives a la contrainte codee (voir section 
2 5 plus bas) 

— soit I'adresse ou le code d'acces d'une zone memoire contenant ef- 
fect ivement les informations relatives a la contrainte codee 

Parmi les informations relatives a la contrainte codee, doit toujours fig- 
urer un code permettant de reconnaitre le type de la contrainte codee 
30 (cf section 1.1.3). 
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• next = 

zone memoire contenant : 

- soit un code indiquant que la contrainte consideree est la derniere 
des contraintes afFectant I'atome A{k). 

- soit I'adresse d'une zone memoire contenant la contrainte suivante 
(voir section suivante). 

Precede de stockage de toutes les contraintes affectant un atome 

Le codage d'un atome A{k) decrit dans la section 1.3.2 contient un champ 
const que nous utiliserons pour stocker : 

• soit un code indiquant que I'atome A(k) n'est affecte par aucune con- 
trainte 

• soit I'adresse d'une zone memoire contenant la premiere contrainte af- 
fectant I'atome A{k) 

Si I'adresse de la n^^^ contrainte relative a I'atome A{k) est stockee dans le 

1 5 champ next de la (n - l)eme contrainte (n > 1), alors, comme le suggere la 

Figure 9 on est assure que le procede de codage decrit permet d'attacher a 
1 atome A{k) toutes les contraintes qui lui sont relatives. Ce mode de codage 
permet de rajouter ou de supprimer facilement des contraintes sans remettre 
en question rorganisation des autres donnees. 

2 0 Exemples d'informations relatives a une contrainte donnee 

Dans la section 1.1.3 nous avons presente a titre d'exemple trois types de con- 
tramtes tres importants pour le codage et I'ajustement des surfaces complexes 
rencontrees en biologie ou geologic. A chacun de ces types de contraintes cor- 
respondent des informations qui doivent etre associees au champ contrainte 
2 5 decrit au paragraphe precedent et schematise sur la Figure 9: 

• Pour une contrainte de type noeud de controle, en plus du type, il est 
necessaire de stocker au minimum les donnees suivantes : 

- les coordonnees ^l, ^J) approchees du noeud <p^ correspondant 
au noyau de I'atome auquel est attachee la contrainte 
30 - le coefficient de certitude (nombre positif) mesurant la confiance 

que Ton peut avoir dans les coordonnees (^f,<^^,y?|). 
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• Pour une contrainte de type lien vectoriel, en plus du type, il est 
necessaire de stocker au minimum les donneies suivantes : 

— Tadresse du noeud (p^ auquel est lie le noeud (p^, correspondant au 
noyau de Tatome A[k) auquel est attachee la contrainte 

5 - les trois composantes approchees du vecteur (cp;^ — ^k) joignant le 

_ noeud (p^ a au noeud 

— le coefficient de certitude (nombre positif) mesurant la confiance 
que Ton peut avoir dans les composantes approximatives du vecteur 

10 • Pour une contrainte de type normale, en plus du type, il est necessaire 
de stocker au minimum les donnees suivantes : 

— les trois composantes approchees du vecteur V orthogonal a la sur- 
face S au noeud (p^ correspondant au noyau de Tatome A{k) auquel 
est attachee la contrainte 

-15 — le coefficient de certitude (nombre positif) mesurant la confiance 

que I'on peut avoir dans les composantes approximatives du vecteur 

V . — 

La presentation detaillee de ces trois types de contraintes est proposee dans 
les sections 2.5 et 2.7.6 de I'annexe, 

20 Variante 

Bien que Ton n'ait parle dans ce qui precede que d'un seul coefficient de 
certitude pour les composantes en x, y et z de chaque contrainte, on pourra 
affecter a ces 3 composantes 3 coefficients de certitude differents si necessaire. 

Codage d'une surface 5 

2 5 Afin d'optimiser la programmation de la variante de la methode DSI presentee 
en annexe, nous proposons de coder une surface S composee de facettes tri- 
angulaires et/ou polygonales sous forme d'un ensemble de N atomes dont 
les noyaux (noeuds) sont les sonmiets de ces triangles et/ou polygones, Les 
donnees et codes correspondant a chacun de ces atomes {A{1)^ . . . , A{N)} 

30 peuvent etre stockes : 

# soit dans un tableau comme le suggere la Figure 10 
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• soit dans le champ atom de zones memoires chainees entre eUes par leur 
champ next comme le suggere la Figure 1 1 .Dans ce ca^, le champ next 
assoae a Tatome aumero k doit contenir I'adresse ou le codHw 
de 1 atome smvant numerote (k + 1); le champ next associe au demie 
aUeTnt ^^ ^^'^^"^"^ atome a"al 

La deuxieme solution est la plus interessante si I'on souhaite pouvoir ajouter 
ou s„er fakement des atomes, par contre elle consomme d'avant^ge d" 
place memoire que la oremiere. mintage ae 



10 



place memoire que la premiere. 

1.3.5 Variantes 



Pour d« commodites de presentation, les zones memoires utiUsees pour coder 

• un nom propre 
15 • une place dans la zone memoire correspondante 

II est bien evident que le procede de codage decrit dans ce brevet ne denenrf 
pa. des noms ni de I'ordre de ces champs a I'interieur desTnes nS^^^^^^ 
chl^'^T'i" ""^"^i ''"^ P°-ble de rajouteTdL 

1.4 Ajustement d'une surface 5 
1.4.1 Introduction 

Etant donne une surface 5 codee suivant le procede decrit dans la section 
1.3 nous nous proposons dans ce qui suit de decrire un procede b^e sur ce 

"efnoel"^!^'""^'-^';"^"''^^ optimale les coordonnees drcitlin 

des noeuds {y>i , - • - , et ceci de sorte que la surface 5 : 

• soit la plus lisse possible au sens du critere de rugosite globale de la 
methodeDSI (voir references [7], [8] et annexe ainsi^ la section l li^ 

• respecte le mieux possible les informations et contraintes connues avec 
30 plus ou moms de certitude concernant la forme de cette surface ( voir 

sections 1.1.1 et 1.1.3) ^ 
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Pour ce faire, on utilisera une variante originale de la methode DSI presentee 
en annexe specialenient con?ue pour utiliser le codage decrit dans la sec- 
tion 1,3. Compte tenu des informations disponibles observees ou mesurees, 
cette methode ajuste successivement et defagon iterative les trois coordonnees 
(S^a> V^a'V^a) chaque noeud Cp^ non fixe suivant une formule du type : 



nouvelle valeur de 
nouveile valeur de 
nouvelle valeur de 



/^( anciennes valeurs de c^^f, 
anciennes valeurs de ^j,. 
/^( anciennes valeurs de (^f,. 



Plus precisement, si I'on se reporte a I'annexe presentant la variante utilisee de 
1 0 la methode DSI, ces fonctions d'ajustement /*(), /|() et Jl{) sont deduites de 
la forme locale des equations DSI au point numero a et ont la forme suivante 
(voir sections 2.4.1, 2.4.2 et 2.7.3 de I'annexe) : 



/^(•••) = 



/«(•••) = 



Ml 



Ml 



{ 

+ E-{rr-<3r}) 



{ 



{ 



E-frr - Qn) 



25 



avec 



' Mr 


= M„+i:7r 






^ Mr 




M„ : 
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La nature et le role des coefficients de ponderation {v-{k)} et u(k) ont ete 
mdiques section 1.1 4 et sont detailies dans I'article [8] ainsi qu 'en annexe Les 
valeurs des termes {Ff" , rr , rf"}, {7^ , 7?°, 7/°} et des termes {Qf^, 
dependent des types et des valeurs des contraintes affectant le noeud i et 
leur formulation exacte est presentee en annexe. 

Dans le cas paxticulier oii les coefficients de ponderation {i;-{fc)> corre- 
spondent a la ponderation haxmonique presentee section 1.1.4 et en annexe 
les fonctions de ponderation /-(), et /^() se simplifient et prennent la 
lorme suivante : 



/.(a) . |A(a)| . ( ^ ^l) 

E /(*)■{( E ^s)- 1 A(i)i-i»d 

fiia) . |A(a)| • ( E V^O 

/(*)■{( E vg)-|A(*)|.rf} 

/.(a) . |A(a)| . ( E Vl) 
E Mfc)-{(E V>'p)-\A(k)\.^l] 

E{rr-gr} ] 



' Afr = M„+E7r 
^:^ = M„-i-E7r 
' ' ^r = A^„+^:7r 

M„ = ^(a).iA(o)|='+ E Mfc) 

V Jfc6A(o) 
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Nous proposons dans ce qui suit un precede de calcul de ces fonctions base 
sur le codage presente au paragraphe 1.3: 

1A.2 Notations 

Afin de simplifier la description des precedes de calculs presentes dans ce qui 
5 suit, il est commode d'introduire les notations et definitions suivantes : 

• Par definition, nous appellerons bloc d'operations, toute iiste d'operations 
precedee d*une accolade ouvrante et terminee par une accolade fer- 
mante Les operations figurant dans un bloc d'operations peuvent 
etre des operations elementaires et/ou d'autres blocs d'operations. 

10 • Par definition, nous appellerons memoire de travail toute memoire utilisee 
pour stocker des result ats intermediaires au cours d'un procede de calcul. 

• Par definition, nous appellerons "'chargement" de I'expression e dans 
une memoire m Toper at ion consist ant a evaluer la valeur de e puis a 
ranger le resultat dans la memoire m. Par convention, une telle operation 

-l 5 sera notee de la fagon suivante : 

m -H- e 

• Pour toute memoire m utilisee pour stocker une valeur numerique, nous 
appellerons "incrementation" de m par Tincrement i Toperation consis- 
tant a calculer la valeur m + i puis a ranger le resultat dans la memoire 

2 0 rri' Pax convention, une telle operation sera notee de la fagon suivante : 

m ^ m + i 
Dans I'operation ci-dessus, Tincrement i peut etre : 

— soit une const ante 

— soit le contenu d'une memoire 

2 5 — soit une expression arithmetique 

Dans Foperation ci-dessus, I'increment i peut etre : 

• soit une constante 

• soit le contenu d'une memoire 

• soit une expression arithmetique 
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1 .4.3 Procede d'a justement dans 1g cas de ponderateurs 
harmoni ques 

Dans le cas oil les coefficients {v°ik)} sont les ponderateurs harmoniques 
presentes dans la section 1.1.4 et dans I'annexe, les valeurs f et r des 
fonctions d'ajustement /-(), et /^() relatives au noeud peuvenretre 
calculees de fagon efficace a I'aide du procede suivant base sur le codage de 



0) Allouer les memoires de travail suivantes 



10 



fx fy f2 

J ay J ay J a 






1) Effectuer les initialisations suivantes (rordre n'a aucune impor- 
tance) : 



= <r5J = = 0 



20 



2) Determiner le nombre de satellites de A{q). Ce nombre pent 
etre obtenu a partir du champ nb-sat de A{a) 

3) Soit A(o;) I'ensemble des satellites de Tatome A(a) auxquels on 
accede par le champ sat de cet atome. Pour chaque satellite ifc € A(a), 
repeter les operations du bloc suivant : 

{ 
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10 



15 



20 



25 



3.1) Determiner le nombre de satellites de'A(k). Ce nom- 
bre peut etre obtenuL a partir du champ nb_sat de -A(ib) 

3.2) Determiner les coordonnees {v^ly^Pk^^D noeud 

de I'atome A{k). Ces coordonnees peuvent etre obtenues a 
partir du champ noeud de A{k). 

3.3) A partir des in&rmations contenues dans A{k)^ calculer 
ou determiner la vafeur du ponderateur fi{k) (voir section 
1.1.4 et annexe). 

3.4) EfTectuer les operations suivantes (I'ordre n'a aucune 
importance), 

<^;x ^^i + 

4 = 5^ = .- = o - 

3.5) Soit A(k) Toisemble des satellites de Tatome A{k) 
auquel on accede par le champ sat de cet atome. Pour 
chaque satellite ^ £ A(t), a effectuer les operations du 
bloc suivant : 

{ 

3.5.2) Determiner les coordonnees actuelles 
i<ppi<p^j^p) du ttoyau tpp de Tatome A(fl). Ces co- 
ordonnees peuvwt etre obtenues a partir du champ 
noeud de A{/3), 

3.5.3) Effectuer ks 3 operations suivantes (I'ordre n'a 
aucune importance) : 



30 



} 



1.2 



3.6} Effectuer les 3 operations suivantes (I'ordre n'a aucune 
importance) : 



n 

fa 



-i? + Mt)x(*^-«itx<pf) 

- fa + M*) X - »fc X <pl) 
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4) EiFectuer les operations suivantes : 



Ma 
Ml' 



10 



15 



20 



5) Soit C(a) Tensemble des contraintes attachees a I'atome A(a) 
auquel on accede par le champ const de cet atome. Pour chaque con- 
tramte c € C(a) repeter les operations du bloc suivant : 

5.1) Determiner le type de la contrainte c 

5.2) En fonction du type de la contrainte c, effectuer les deux 
operations suivantes : 

5.2.1) Extraire les informations specifiques de la con- 
trainte c 

5.2.2) En fonction des informations extraites en (5.2.1), 
calculer les valeurs des expressions suivantes (voir sec- 
tions 2.5 et 2.7 de I'annexe) : 

pxor -pyot -nza 

5.3) Effectuer les operations suivantes (Pordre n'a aucune 
importance) : 



25 



Ml' 

M^y 



/« + (rr-Qf«) 

fi + (rr° - Qr ) 

/a + (rr-Qf«) 



30 



6) Effectuer les operations suivantes (I'ordre n'a aucune importance) 



J a 
J a 



(p(of) X n„ X o-J - fS)/Ml' 

(/*(af) X Ha X o-v - fy)/Miy 

ift(a) X n„ X (T^ - f^)/M, 



A la fm de ce procede de calcul, les memoires de travail /v et /| contiennent 
les valeurs des fonctions d'ajustement /^(), fvQ et flQ utilisees pour ajuster 
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les valeurs des coordonnees du noeud tp^ : 




5 Le precede d'ajustement decrit ci-dessus ne doit etre applique qu'aux noeuds 
dont la position est non fixee, c'est-a-dire aux noeuds qui ne sont pas des 
noeuds de controle (voir section 1.1.3). 



1*4 .4 Vaxi antes 



Variante 1 



1 0 La mise en oeuvre du procede d'ajustement (base sur la methode DSI) presentee 
au paragraphe precedent dans le cadre de ponderateurs harmoniques peut 
facilement etre adaptee dans le cas de ponderateurs quelconques. La structure 
globale du procede reste la meme, seule change ['incrementation des difFerentes 
memoires de travail allouees a I'etape (0); la fagon d 'increment er ces memoires 

1 5 de travail depend des ponderateurs choisis. 

Par exemple, si nous choisissons les ponderateurs '{v°{k)y suivants, 




si a = k 
si a € A(ft) 



FEUILLE DE ^EMPLACEMENT 



3NSDOCID: <WO. 



.9104544A1 J_> 



wo 91/04544 



- 28 - 



PCr/FR90/00676 



/»(•••) = 



1 0 



alors on verifie facilement que les fonctions d'ajustement associees sont : 

[ /'(^)-( E ^f)/|A(a)| 

[ E 'Pl)/\Hc.)\ 

- E{rr-Qr} J 

^fc€A(a) ^ 

E{rr-Qn ] 



15 



avec : < 



f ^r = M,+E7; 



,za 



Ma = /x(q) + {/^ W/|A(ifc)p} 
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0) Allouer les memoires de travail suivantes 



fy p 



1) Effectuer les initialisations suivantes (I'ordre n'a aucune impor- 
tance) : 



10 



15 



2) Determiner le nombre tia de satellites de -A(af). Ce nombre pent 
etre obtenu a partir du champ nb^at de -A(Qf) 

3) Soit A(of) Tensemble des satellites de Tatome A{a) auxquels on 
accede par le champ sat de cet atome. Pour chaque satellite k £ A(Qf), 
repeter les operations du bloc suivant : 



{ 
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3.1) Determiner le nombre de satellites de 'A{k). Ce som- 
bre peut etre obtenu a partir du champ nb^at de A(k) 
Jr/ffectuer ensuite I'operation suivante : 



1 0 



15 



20 



25 



30 



35 



3 2) Deterrr^ner les coordonnees (^%,<pl^l) du noeud <p, 
de 1 atome A(*). Ces coordonnees peuvent etre obtenues k 
partir du champ noend de A{k). 

3.3) A partir des informations contenues dans ^(Jb), calculer 
ou determiner la valeur du ponderateur f,{k) (voir section 
1.1.4 et annexe). 

3.4) Effectuer les operations suivantes (I'ordre n'a aucune 
importance). 



»5 



3.5) Soit A(*) Tensemble des satellites de I'atome A(k) 
auquel on accede par le champ sat de cet atome. Pour 
chaque satellite /3 € A(k), /3 ^ a effectuer les operations du 
bloc suivant : 

{ 

3.5.2) Determiner les coordonnees actuelles 
i'P0,<p0,<p'p) du noyau ip^ de I'atome A{/3). Ces co- 
ordonnees peuvent etre obtenues a partir du champ 
noeud de A(^). 

3.5.3) Effectuer les 3 operations suivantes (I'ordre n'a 
aucune importance) : 



} 



3.6) Effectuer les 3 operations suivantes (I'ordre n'a aucune 
importance) : 



J a 

n 



fl + m(A) X {,yjnl - ^l/n^) 
r. + Kk)x(s'^/nl-.<pl/n,) 



FEUILLE DE REMPLACEMENT 



wo 91/04544 



- 31 - 



PCr/FR90/00676 



} 



4) EfTectuer les operations suivantes : 









Ml' 




Ma 






Ma 






Ma 



10 



15 



20 



25 



30 



5) Soit C7(ar) I'ensemble des contraintes attachees a I'atome ^(or) 
auquel on accede pat le champ const de cet atome. Pour chaque con- 
trainte c € 0{<x) repeter les operations du bloc suivant : 

{ 

5.1) Determiner le type de la contrainte c 

5-2) En fonction du type de la contrainte c, efTectuer les deux 

operations suivantes : 

5.2.1) Extraire les informations specifiques de la con- 
trainte c 

5.2.2) En fonction des informations extraites en (5.2.1), 
calculer les valeurs des expressions suivantes_(voir sec- 
tions 2.5 et 2.7 de Tannexe) : 

pxa pyc pra 

5.3) EfTectuer les operations suivantes (Pordre n'a aucune 
importance) : 



} 



Mr 



Mr 



Ml' ^ Ml' + Tf " 



6) Eifectuer les operations suivantes (Fordre n*a aucune importance) 



n 
fx 
n 



X ^Jn, - m/MT 
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Les parametres numeriques stockes dans le champ info des atomes peuvent 
egalement etre interpoles par la methode DSI (yoir section 2.1.2 et 2 2 de 
I annexe). Dans ce cas, la structure du procede d 'interpolation est quasi- 
ment identique a celle du procede d'ajustement de la surface presente ci- 
dessus; le seul changement consists a ne conserver qu'une des trois fonction. 
d ajustement normalement dediees aux trois coordonnees des noeuds et a cai- 
culer celle-ci en fonction du parametre interpole. 

La pr6sente invention trouve application tout 
particuli^rement dans les domaines suivants : 

- En g^ologie et en g6ophysique, elle permet de 
crier un module de surface, par exemple de construire une 
maquette solide en trois dimensions, sur laquelle on peut 
effectuer une simulation g6ophysique ou une simulation de 
1 -exploitation de richesses souterraines ; elle permet 
15 Egalement d'obtenir un ensemble de donn6es propres A 
initialiser un simulateur num6rique ou analogique de 
ph6nom^nes giophysiques ou d • exploitation de richesses; en 
outre, elle permet, h 1 • aide d'une station de tralail 
graphique, d'afficher sur un 6cran ou d'imprimer sur un 
20 support des images permettant de contrdler les simulations 
ci-des8U8 ; 

-en biologie ou en m^decine, 1' invention permet de 
construire des modules en trois dimensions d'organes et/ou 
de prothdses, ou encore de simuler les effets d' operations 
25 en chirurgie plastique; de m%me, des images en deux 
dimensions affich6es ou imprim6es permet tent ici encore de 
contrdler les operations ci-dessus. 

Mais il est clair que 1' invention s' applique plus 
g6n6ralement ^ la mod^lisation de surfaces de type 
quelconque, et notamment de toute surface naturelle ou 
industrielle. Elle peut trouver un int6rit tout 

particulier dans le cadre de la conception assistfee par 
ordinateur (CAO) . 

Enfin, la presente invention n'est nullement 
liBit6e 4 la forme de realisation d^crite ci-dessus et 
representee sur les dessins, nais I'homme de I'art saura y 
apporter toute variante ou modification conforroe k son 
esprit . __, 
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ANNEXE 



2. NOUVELLE VERSION DE LA METHODE DSI ADAPTEE A LA 

MODELISATION INTERACTIVE DES SURFACES COMPLEXES. 

Dans les domaines de la Geologic et de la Biologit, il est frequent d 'avoir a modeliser des 
surfaces complexes correspondant par exemple a des interfaces entre zones de nature 

5 ou de proprietes differentes. Les techniques classiques de modelisation basees sur les 
interpolations de Bezier ou sur les fonctions Splines (voir ref. [9]) sont tres mal adaptees 
a la prise en compte de ces donnees beterogenes* Nous proposons ici une approche 
differente basee sur la methode DSI (Discrete Smooth Interpolation - voir ref, [8]), Dans 
cette approche t les surfaces sont modelisees avec des facettes triangulaires irregulieres 

10 dont les sommets doivent etre determines afia de tenir compte du plus grand nombre 
de donnees heterogenes possible. 



15 2.1 Introduction : 

Presentation des problemes etudies 

2.1«1 Surface S definitions associees 

20 Comme le montre la Figure 1 , soit S une surface composee de facettes triangulaires ad- 
jacentes plongee dans I'espace euclidien 3D (0|x, y, z), Cette surface n est pas forcement 



9104544An I > 
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10 



n'est pa^ forcement connexe ni fermee. Dans ce qui suit, nous utiliserons les 
notations suivantes : 

G = G{S) = ensemble des cotes des triangles 
^ n = ^{S) = ensemble des sommets des triangles 
, AT = |n| = nombre de points de 

G est un graphe dont les noeuds sont les points de et nous identifierons cet 
ensemble avec ensemble des indices utilises pour les numeroter. 

Le voxsmage N{k) d'un noeud Ar € Q est defini par : 

N{k) = sous ensemble de Q constitue par les noeuds de 
n qui peuvent etre atteints a partir du noeud k 
en s{k) pas au plus ou s{k) est un nombre positif 
donne. 

^^mJ^ """i!' P°^^^^PP°s«^°°s ces voisiaages N^k) possedent la pro- 
priete de symetrie suivante : ^ ' f 

Les voisinages ainsi definis engendrent une topoiogie sur G qui peut etre con- 
sideree comme une approximation de la topoiogie de S elle meme. 

2.1.2 Profeleme 1 : Estimation d'une fonction <p definie 
sur Q 

20 Supposons que la position de chaque noeud Ar € est connue et soit <p(k) 
une fonction definie pour chaque noeud A: € n et connue seulement sm un 
sous-ensemble de n : 

L = ensemble des noeuds I e n oil ip{l) = est connu 
/ = ensemble des noeuds i^n ou ^{i) = est inconn 



25 



_ ^ ^ rv , V'l est inconnu 



Generalement, L est different de 0 et 11 y a une infinite de fonctions ^(k) 
defimes sur fi et interpolant les valeurs {^(l) = ^, : / g j^j. N^^re objectif est 
de selectionner parmi toutes ces fonctions celle qui minimise un critere donne 
■nr{<p) mesurajit : 

30 •la rugosite globale de la fonction (p(-) 

• I'ecaxt entre la fonction ^{.) et un certain nombre de contraintes et 
donnees qu'elle doit satisfaire 
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2.1.3 Probleme 2 : Ajustement de la surface S 

Dans ce cas, seuls les points d'un sous ensemble L de noeuds A: 6 fi ont une 
position connue ip{k) = (pj^ : 

L = ensemble de noeuds / 6 ou = (pi est connu 

/ = ensemble des noeuds i ^Vt ou ^{i) = <p- est inconnu 

== n-L 

Comme dans le cas du premier probleme, L est generalement different de Q 
et il y a une infinite de surfaces S interpolant les points '{<p{l) = : I E L"}. 
Notre objectif est de selectionner parmi toutes ces surfaces, celle qui minimise 
1 0 tin critere donne mesurant : 



# la rugosite globale de la surface S 

# Tecart entre la surface S et un certain nombre de contraintes et donnees 
qu'elle doit sat isf aire 

Compte tenu de I'independance des composantes (ip'{k),ip^{k), ip^(k)) de chaque 
15 vecteur <p{k)y il est facile d'utiliser la solution du premier probleme pour 
resoudre le second en posant simplement : 



^ — Dans le cas tres particulier ou Tensemble Q correspond aux noeuds d'une 
grille rectangulaire posee sur la surface, on peut songer a utiliser une methode 
20 d 'interpolation du t3rpe spline, Briggs ou Bezier (Voir ref. [9], [1] et [3]); 
malheureusement ces methodes sont inappliquables dans le cas des grilles 
irregulieres et ne permettent pas de tenir compte des contraintes presentees 
aux paragraphes 2.2.3, 2.5 et 2.7.2. 



2*2 Interpolation 

d'une fonction (p{k) definie sur Q 

Dans cette section, nous introduisons brievement la methode DSI (Discrete 
Smooth Interpolation) presentee dans Tarticle [8]. 
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15 



20 



25 



30 



2.2.1 Notations ~ 

Dans ce qui suit, nous noterons ^ une matrice colonne de taiile N telle que : 

<Pk = VA; 6 n 

Notre probleme ne depend pas de la methode de numerotation des noeuds de 
la gnlle Afin de simpHfier les notations, nous supposerons queTes^nt 
de la matnce ^ ont ete eventueUement permutes de telle fa,on que i Tuisse 
etre partagee en deux sous-matrices <pj et telles que . ^ "^"^ P"'^^^ 



= lK] 



V/ = 



'PL - 



avec 



€ / Va 



L J 



avec lp£ L ^ 



De plus nous noterons I| • ||^ une semi-norme associee a une matrice carre. 
^ma-defime positive [D] de taiile (^V x N) de telle fa^n qu^pTur to^^ 
matnce colonne X de taiUe on ait : m , pour louce 

2.2.2 Definition d'un ciitere de rugosite globale R{^) 
•[r (fc)} sont des coefficients de ponderation positifs, negatifs ou nuls donnes : 

or6W(it) 

i^Cv^lA:) pent etre utiUse pour construire un critere de rugosite globale Ri^^ 
^^. t^^^^^"" '"^^"^'^ °" '^'^ ponderaLn doniri 

k 

i2(<^) est completement defini pax la donneedes coefficients {t;«(fe)} et A^fA:) et 
il est facile de venfier (voir ref. [8J) qu'il existe toujours une matr ce x ^ 
symetnque semi-definie positive [W] telle que : V ^ -'v ; 

Ri<p\k) = <p* . [W{k)] . <p 
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2.2.3 Prise en compte de contraintes lineaires an sens 
des moindres carres 

Considerons une matrice carree [A^] de taille {N X A^) et une matrice colonne 
Bi de taille A'' definissant la contrainte lineaire suivante : 

Etant donnee une matrice semi-definie positive [Z?^], le symbole figurant 
dans la formule ci-dessus a la signification suivante : 

II [A^ • <p — Bi ||£). aussi petit que possible 

S'il existe plusieurs conditions de ce type, nous proposons de mesurer le degre 
10 de violation de ces contraintes par le critere p{^) defini comme suit : 



15 



20 



2.2.4 Solution du probleme _ 

Parmi toutes les fonctions^(A:) definies sur CI et interpolant les donnees 
{(p{l) = ipi : I £ L), nous en choisirons une qui minimise le critere suivant : 

En developpant i?*(v)> on obtient : 

i?*(v?) = <f''[W*]'tf> - 2 . <i?*.[Q] + C 

Q = E W ■ W ' 

i 

c =i:iiB,iii, 

La partition de la matrice 9 induit une partition similaire des matrices \W*] 
et Q utilisees pour definir ' 



avec < 



"Pi 
9l 



LL 
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L^coadition aB-(s=)/av=, = [0] implique I'equation suivante apoel« «^„».- 

2.3 Unicite de la solution de I'equation DSI 
Definition 

appartenant k L. ^ contient au moins un noeud 

theoreme 

Si le critere de rugosite globale R(^) est tel que : 



4) t;*(ifc) < _ ^0 ^ 

alors. I'equation DSI basee sur R^^) admet une solution unique. 

Demonstration 

20 Le theoreme ci-dessus ne presente an 'itn inf^-Af • 

cidence sur k precede dfcit d^ k^r^^t btettrc W ~ 
doimerons pa, ici sa demonstration. ' ^' °<"^ 

2.4 rorme locale de I'equation DSI 

da^ns le presLat "«-"««t precisement ceUe qui est utilisee 
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Afin de simplifier les notations, nous supposerons dans ce qui suit que les 
matrices semi-definies positives [Z)J intervenant dans Tequation DSI sont 
diagonales. De plus, lorsqu'il y a des contraintes lineaires du type [A-] -(^ ~ 
nous utiliserons les notations suivantes : 

{licne k 
colonne h 

Qf = element de la matrice colonne ([A-]*[i),]B.) 
2.4.1 Calcul de SI^{<p)/d<Pa 

Compte tenu de la definition de R{ip\k), on peut ecrire Timplication suivante : 
10 = E + E E ^^ik)v'r{k).^,,p^ 



0^1 



' 2 . K(ifc))2 • <,3„ + 2 . . E ""^^k) • si a € N{k) 



, 0 sinon 



On deduit de cette relation que : 

1 ail(y) 1^ dRim _ 1 ^ dR{ip\k) 

= Va- E IJ^{k){v%k)f + E {Kk)v''{k). E ^^W-V^ 

Compte tenu de la structure diagonale de [i?,], il est facile de verifier que le 
^eme ^i^nient de a || [A,] • - -B,- jd^^ est tel que : 

20 = -E A^' (A^")' + E{a-"A^*= • E Af'^v^ - Q° 

k k ^ ^ 
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Compte tenu de la definition de if'Cv), on en deduit que : 



1 R^jy.) 
2' d<p, - 



- EQr 



La solution y,, est teUe que dlf{v)/d<p, = 0, done la composante :^ 

de (f doit satisfaire I'equation DSI(a) suivante : 



10 



DSI(a) : 



+ E{rr-Qf}} 



avec : ^ 



Tf = EA-*(Art^ 
<3? = E^f°A-*5? 



15 Nous proposons d'appeler cette equation forme locale de Vequation DSI au 
noeud a. 

2.4.2 Proposition pour un algorithme iteratif 

La forme locale de Tequation DSI presentee dans la section precedente suggere 

pour estimer la solution Par exemple, a Titeration 
i?ct/^^\^^ composante 9L"-^'> de la solution ^i^^^) doit satisfaire I'equation 
DSI(a) ce qui permet de proposer le procede iteratif suivant : 

soit / I'ensemble des indices de noeuds ou <p^ est inconnu 
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soit ip une solution approchee initiate 

tant.que( il est necessaire de faire de nouvelles iterations ) 

{ 

poiix.to"at( ct € / ) 

{ 

= {.••} 

} 

Comme on peut le voir, cet algorithtne tres simple n'utilise pas exp licit ement la 
matrice [W*] intervenant dans Tequation DSI mais il requiert de recalculer un 
grand nombre defois des produits du type {t;°'(fc) •t7^(fc)} qui sont precisement 
ceux utilises pour construire [IV]. 

1 0 En fait, si la solution approchee initiale est proche de la solution exacte, 

alors il n'est pas necessaire d'effectuer beau coup d 'iterations et le surcout de 
calcul devient negligeable. Par exemple, cette situation se presente dans un 
contexte interactif ou la solution initiale est prise egale a la solution avant 
que quelques modifications locales soient efFectuees pax I'utilisateur. Dans ce 

15 cas, malgre le leger surcout engendre, la fonne locale de Tequation DSI est 
preferable car elle est bien plus facile a mettre en oeuvre que la fornae globale. 

Afin de montrer que I'algorithme propose converge, on peut remarquer 
que R'isp) peut etre exprime sous forme d'une fonction de 9^ ayant la forme 
suivante : 

— ir(v.) = A.<^^ + J9.<^, + c 

Les coefficients jB et C sont independants de ip^ et on a : 

Le minimum du critere global R*'{(p) est atteint pour la valeur ip^ = -B/{2A) 
25 qui correspond precisement a la valeur fournie par I'equation DSI(a). Ceci 
nous permet de conclure que notre algorithme converge car, a chaque iteration 
du processus iteratif, la valeur de la fonction positive ou nuUe R^{^) decroit : 
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2.5 Prise en compte de donnees imprecises 

Dans ce qui suit, nous introduisons deux types de contraintes Uneaires qui sont 
t:ZZ7Z T^T^^:- 7^^"-*-° g-metrique ainsi qu'on I verra 
5 donnees imprecises concernant les valeurs inconnues {^x : X e 1} De telles 
donnees sont supposees associees a des coefficients de certitude ^ R+ qui 
sont en fait des coefficients de ponderation positifs. 

2.5.1 Cas d'lme donnee ponctueUe imprecise 

. n ^""'^ ^ % f f'^'^'f ''''^'''^ ^^^"^^ 1^ ^^^^ inconnue est supposee 
1 0 prodie dW donnee incertaine ^,avec un coefficient de certitude donnlegS 



'A 



ou ~ correspond a la condition suivante : tp^ : 

^\ ' \ fx — <Px\^ minimum 

15 Une telle con^ion est aisement prise en compte par la methode DSI sous 
torme d une contrainte avec : 



If = / +1 sik = k = \ n* _ f <^A si 
[ 0 smon ~ \ 0 si 



k = X 
sinon 



et £)*^ = i si k = h = X 
20 ' 0 sinon 



DSir^^^'^r? ^h^'- <^o-respondants utiUses dans I'equation locale 

Dbl(a) sont alors definis par : 



Tf = 0 Qo ^ f ' <Px si a = A 
' \ 0 sinon 



25 et 7? - J si a = A 



,2 

Sinon 
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2.5*2 Cas d'line donnee differenti^lle imprecise 

Soit A € / un index de noeud pour lequel la valeur inconnue (p\ est supposee 
liee a une autre valeur correspondant a un autre index ^ X. Supposons 
que ce lien ait la forme suivante ou A;^^ est une valeur donnee : 

5 i^^-^x)^ 

Cette relation est supposee imprecise et est interpretee de la fa^on suivante 
ou est un coefficient de certitude : 

°l^'\{^l-<Px)- ^xj^ minimum 

Une telle condition peut etre facilement prise en compte par la methode DSI 
10 en tant que contrainte avec : 



15 



' -1 


si /: = 


h = X 




-1 


si A: = 


h = ft 




< +1 


si /: = 




= A* 


+1 


si fc = 


PL et h 


= A 


0 


sinon 















+ A 



SI 



- s\k = 

0 sinon 



et D^^ = I i 



0 



si A: = A = A 
si A: = A = 



smon 



Les coefficients Ff , Qf et 7f correspondants utilises dans Tequation locale 
20 DSI(a) sont alors definis par : 



rr 



-s'am "pI si q = a 

-S'L • si or = /X 

0 sinon 



si a = A 
0 sinon 



25 



et 7? = 



,2 



t37 



Am 



si a = A 
si a = ^ 
sinon 
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2.5.3 Choix des coefficients de certitude 

Considerons le terme de I'equation DSI(a) : 

avec : f,(k) {v''{k))^ 

Dans les deux exemples presentes ci-dessus, les termes -v? sont .nif « 
zero, soit correspondent au coefficient de certitude ' ^^^""^ ^ 

7f = 

1 0 M pourcentage donne p, de 

erf = • Af„ Pi > 0 
Le terme M* de I'equation DSI(a) est alors egal a : 

= M„ . (1 + + . . .) 

oil p.(a) est soit egal a un pourcentage donne p. ou egal a zero. 

15 2.6 Choix des coefficients de ponderation 

Le choix des coefficients de ponderation {t;«(jb)> et 4u(m est r^rr.^rf 
hbre excepte que les coefficients ^k)} IJeJAlTo^^ 

^TicTeZ ^^'^^ ^^^"^^^^ ^^^^ '^ - deu7fa™de 

20 2.6.1 Choix de coefficients {v°'{k)} : 
ponderation harmonique 

Soit A(A:) T'orbite" de k definie de la fagon suivante : 

A(k) = N{k) ~ {k} 
Soit |A(A:)( le nombre d'elements de A(jb). Nous dirons aue no,,. „tn^c 

r^rA:) = / si a = k 

\ 1 si a 6 A(A:) 
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Les fonctions harmoniques ont la propriete caracteristique d'etre egales en tout 
point a leur propre moyenne sur un cercle centre en ce point. Comme on pent 
le voir, R{<p\k) = 0 si (^j^ est egal a la moyenne des valeurs de entourant 
le noeud c'est pourquoi, par analogic avec les fonctions harmoniques, nous 
proposons d'appeller harmoniques les coefficients {v^Xk)} presentes dans cette 
section. 



10 



Si nous developpons I'equation DSI(a) avec ces coefficients, nous obtenons 
la formule suivante qui pent etre facilement transcrite dans un langage de 
prograrrmiation : 



DSI(a) : 



jfe€A(a) fl€A(fc) 



avec 



fc6A(a) i 



15 2.6.2 Choix des coefficients {fi{k)} 

Les coefficients {ft{k)} ont ete introduits afin de moduler localement le lis- 
sage de Tinterpolation. En rabsence de raison particuliere, nous suggerons 
d'utiliser une 
ponderation uniforme : 

/i(fc) = l vfcen 



20 



25 



Panni toutes les ponderations non unifonnes, la ponderation suivante semble 
particulierement interessante s'il est necessaire d'avoir une fonction lisse au 
voisinage des donnees : 

''^^^ - \m>lsi keL 
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2.7 Application a la 

Modelisation Geometrique 

meinoae Ubl atm d estimer la surface trianculee S ell^ Tr.«rr,« « i ,, 
: on suppose disposer des inforxnations partieU^t if a^f^ '^^"""^ 

• la position exacte de quelques sommets de triangles, 

• la position approchee de quelques sommets de triangles, 

• des contraintes vectorielles precisant la forme de S, 

• etc . . . 

0 Pour cela on est amene a definir ^„ comme le vecteur courant joi^nant 
1 ongxne de R3 au noeud courant de f2 et nous noterons v^^^et J" sL 
composantessurlabaseorthonormee{ox,oy,or}deR3: ^ 

2.7.1 Defmition d'rni critere de rugosite globale R(<^) 

5 Suivant notre habitude, I'ensemble des noeuds de la griUe est oartacre «n 
deux sous-ensembles / et Z partage en 

{<Pi est inconnu Vie/ 
<?i est connu I ^ L 

Les noeuds / € L dont la position e est connue seront appeles noeud 
tl 7 f f ' °" r determiner la position des poL rrest!att 

M T ^"^J" "'^'^^^ ^^^^ ^(^*) d^-i ci-dessous et bie' st un 

ensemble donne de ponderateurs {V'ik)} soit aussi petit que possible 

^(^*) =11 E ^''W'f^air 

a€Ar(ik) 

Ce critere constitue la forme vectorieUe du critere DSI local. Le critere DSI 
5 global assooe est aJors defini de la fa,on suivante ou les A:)} sont dL 

coefficients de ponderation non negatifs donnes : 

-k 
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Si I'on definit i2((^*), Ri^p^) et R{(p') ainsi qu'il a efce fait dans les paragraphes 
precedents : 

alors, en utilisant le theoreme de Pythagore, on verifie facilement que Ton a : 

2.7.2 Prise en compte de contraintes lineaires au sens 
des moindres c aires 

1 0 Considerons les matrices suivantes dimensionnees de telle f agon que represente 
le nombre total de noeuds de Tensemble CI : 

[D^J . est une matrice diagonale non negative de taille N x N 
\i V — x^y^z : \ [A^J\ est une matrice carree de taille N X N 

BY est une matrice coloime de taille 
•1/ 

15 Si nous voulons que les composantes if^^ ip^ et if^ satisfassent les contraintes 
suivantes 

[AfJ-s^- ~ Bl 
2 0 alors on est amene a introduire le critere sui vant 

f P{.^') = E ll[>ifJ-v'-5f,llLr.. 

Pi^n = Z \\[AU'^'-Bi\\i^, 



avec 
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2.7.3 Solution du probleme 

Parmi toutes les surfaces S interpolant les donnees {3(1) ^ <3, • I ^ L\ 
nous proposons de selectionner celle qui minimise le critL suivaLf : ^' " 



10 



avec : < 



15 



Compte ten. de I'independence de. con>pos»tes ^. et oa en condut 

B^{<p) minimum 



/?*(^*) minimum 
minimum 
, ^{}p') minimimi 



2.7.4 Remarque 

won unique des lors qu il y a au moms un sommet de facette av«r,f 

position connue sur chaque morceau. "^^^ 

2.7.5 UtiKsation interactive 

" s^Lndf ''f S^metrique, I'objectif est de definir interactivement la po 
sition dans 1 espace du graphe G associe a I'ensemble des noeuds n Pour^ 

"::ret=:rnr^^^^^ 

de controles et donnees impTrses^ Zd fi "^^^^^^es (noeuds 

30 aux noeuds a e / en Ta mod rfications sont alors propagees 

* t J en uiiiisant la forme vectorielle de la methode DSI. 
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Nous avons montre dans les paragraphes precedents que ce probleme se 
decomposait en trois sous-problemes correspondant aux trois coordonnees 
{x,y,z) et la meilleure fagon de resoudre ces trois problemes est d'utiliser 
Talgorithme iteratif associe a la forme locale des equations DSI correspon- 
5 dantes- 

D'une etape a Tautre du processus de modelisation on peut modifier non 
seulement la position des coordonnees des noeuds de controle mais Tensemble 
L lui meme peut etre change si necessaire; par exemple, si nous voulons appli- 
quer DSI seulement a un sous ensemble T de noeuds, alors nous devons choisir 
10 L de telle fagon que i D F ou F designe Tensemble complementaire de T dans 

a. 

De plus, compte tenu de la decomposition de notre probleme en trois sous- 
problemes resolus independemment et correspondant aux trois coordonnees, 
on peut, si necessaire, defi.nir des sous ensembles difFerents L^, Ly et Lz pour 
15 chacune de ces coordonnees; par exemple, si les coordonnees (x,y) ne doivent 
pas etre modifiees, alors il suffit de prendre et Ly egaux a I'ensemble Q 
tout entier. 

' 2.7.6 -Prise en compte de doimees geometriques imprecises 

Initialement, la methode DSI a ete congue pour modeliser les surfaces com- 
20 plexes rencontrees dans le domaine des sciences naturelles comme par exemple 
la geologic et la biologie. Dans ce cas, on doit generalement prendre en compte 
de nombreuses donnees impresices comme par exemple : 

( m des noeuds de controle flous 
I • des donnees vectorielles floues 

25 Noeuds de controle flous 

Par definition, nous appellerons "noeud de controle flou" tout noeud A € 
/ dont la position tp^ est connue conune etant approximativement egale a 
un vecteur donne et ceci avec un certain degre de certitude s?^. En se 
referant au paragraphe 2.5.1, on voit que Ton peut prendre en compte une 
30 telle information imprecise a Taide de la methode DSI en introduisant les 
contraintes suivantes dans les criteres i?*(v?*)» ^(v?*) • 
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Contraintes vectorielles floues 

Dans de nombreuses applications, il est necessaire de controler 1« v... 

(v.-9A)joignant deux noeuds^, etc? d'unesurfax:e pJ^^fi ^ T 
- contrainte sera appellee W Jnte ve^or^Uell^^^^^^ T^-^'' 
5 avec un certain degre de certitude n ?! , ^^^^ satisfaite 

ce type sent presentes surTa n^,e ^'^^"^^^ 

En pratique, au moins un des deux points <p. on <? a une nosition ir, 

par exemple celui correspondant a A € / et narmflp. ^ ^^connue, 

citer les deux suivants : ^ ^ ^ ^ > ^t parmi les cas importants, on pent 

10 • Le premier cas correspond aux Fieures 5 » At o k i» j ' • 
la forme d-.n objet de teUe f^^Z ^ 1 SjTo t iflT 
donne A " '^'^ ^ "° vecteur 

a un objet de teUe fagon que chaque _ ,^it orthogonal a un 

= '''' coefficient de ponderation egal a 

Pour le second cas, le probleme n'est pas tres different car on a : 

Dans ce cas nous suggerons de prendre les valeurs des termes (cpf 
(^.^<Pl,^l) mtervenant dans les mexnbres de droite de ces equtd'o'^slgilL 
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aux valeurs correspondantes obtenues a Titeration precedente; pour se ramener 
au cas precedent, il suflBt alors de poser : 

Remarque 

Dans la prise en compte des donnees imprecises presentee ci-dessus, nous 
n'avons utilise qu'un seul coefficient de certitude w^^ commun aux trois com- 
posantes de ces contraintes. En pratique, ceci n'est pas une obligation et, 
T 0 compte tenu de I'independance des trois composantes, il est toujours possible 
d'utiliser trois coefficients de certitude difFerents tn^J, ct^^ et ct^J^ correspon- 
dant a ces trois composantes. 

2*8 Generalisations 

Avant de conclure, mentionnons quelques generalisations evidentes : 

15 •La methode DSI peut etre utilisee avec n'importe quels types de facettes 
polygonales. En pratique, nous conseillons d'utiliser des triangles car ce 
sont les p_olygones les plus simples et les plus faciles a manipuler dans 
un ordinateur. 

— •La methode DSI peut etre utilisee pour estimer des courbes polygonales 
2 0 dans un espace 3D; dans ce cas, les triangles sont remplaces par des 

segments. 

2-9 COlNTCLUSIOlSr 

La forme locale de Tequation DSI constitue un outil simple et puissant pour 
modeliser les surfaces complexes que Ton peut rencontrer par exemple en 
25 geologie ou en biologie. En particulier : 

• Le maillage utilise pour modeliser une surface peut etre completement 
libre et il est possible d'utiliser des algorithmes automatiques pour les 
construire; par exemple,^i le naaillage est compose de triangles, on peut 
utiliser des algorithmes derives de la methode de Delaunay (voir ref. 
30 [2] et [5]). De plus, la taille du maillage peut facilement etre ajustee 

localement en fonction de la complexite de la surface. 
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L r ' ' P"^"**^^ ^ grand nombre 

de donnees precises ou imprecises. Cette possibilite est paTticuiierement 
mteressante dans le doxnaine des sciences naturelles ou^ I'obj J^™ 
pa^ de generer des surfaces esthetiques mais d'ajnster ces surfal a^x 
donnees precises et imprecises de fa^n coherente et efficace. 

• L'algorithme derive de I'equation DSI locale est suffisamment raoide 
pour permettre une utiUsation interactive. ^ 

Nous n'avons rien dit sur la representation des facettes elles-memes Selon les 
appbcations, on pent prendre des facettes planes mais, si cela est Sees" 
est toujours possibk de representer ces facettes pax d;s morceaux de su^^^^ 
non planes mterpolant les noeuds du maillage (voir ref [4]) 
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REVEMDICATIONS 

1. Proc6d§ d'obtention d'une mod^lisation d'une 
surface (S) representative notamment de 1' interface entre 
deux zones de natures ou de propri6t6s diff6rentes dans 
un corps ^ trois dimensions tel qu'une formation 
g§ologique ou un corps vivant, caract6ris6 en ce qu'il 
comprend les 6tapes consistant k: 

- obtenir au moyen d'appareils de mesure un 
ensemble de donn6es g6om6triques relatives A la surface 
et associ6es respect ivement A des points donnas de ladite 
surface? 

- effectuer un maillage de la surface de telle 
sorte que 1' ensemble desdits points donnas soit un sous- 
ensemble des noeuds {<^^ du maillage; 

- m6moriser une adresse donn6e d'une m6moire, 
sp6cifique ^ chaque noeud ( k ^ maillage, les 
informations suivantes: 

* les coordonn6es (Cpj^, ^fc, du noeud consid6r6r 
^ * le nombre (nb-sat) de noeuds satellites du 

noeud consid^r^r 

* une information (sat) permettant d'acc§der aux 
adresses sp6cifiques desdits noeuds satellites et par 
suite aux informations relatives k ces derniers, 

* le cas 6ch6ant des donn6es g6om6triques 
(const) asfioci^es audit noeud consid6r6, 

- pour chaque noeud du maillage r d6finir un indice 
25 de rugosit* local (R(Cf|k)) obtenu k partir d'une somme 

pond6r6e des coordonn6es actuelles du noeud et de ses 
satellites, 

- d6finir la somme (R* «^) ) d'un indice de rugosit6 
global (R(Cp)), constitu6 par 1* addition des indices de 

30 rugosit6 locaux associ6s k chaqfue noeud, et d'un indice 
global (^{^) de violation desdites donn6es g6om6triques ; 

- ajuster de fagon iterative les coordonn6es de 
chacun des noeuds dont on ne connalt pas les coordonn6es 
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precises, en utilisant 4 chaque fois 1' addition d'une 
combinaison pond6r4e des coordonn6.es actuelles des 
satellites et des satellites des satellites dudit noeud 
et d'une combinaison des donn6es g6om6triques associ6es 
5 audit noeud, de mani^re k amener ladite somme k un 
minimum; 

- cr6er une repr6sentation de la surface h partir 
des coordonn6es ajust6es (Cp*C|»acp,5 de chacun des noeuds. 

2. Proc6d6 selon la revendication l, caract6ris6 
10 en ce qu'on effectue un maillage triangulaire de la 

surface. 

3. Proc6d6 selon I'une des revendications 1 et 2, 
caract6ris6 en ce que l'6tape d'ajustement comprend trois 
sous-6tapes distinctes mettant en jeu respect ivement les 

15 trois coordonn6es spatiales Cp,?,q>f) des noeuds et des 

satellites. 

4. Proc6d6 selon I'une des revendications 1^3, 
caract§ris6 en ce que les donn6es g6om6triques 
comprennent des donn6es de coordonn6es des noeuds. 

5. Proc6d6 selon I'une des revendications 1^4, 
caract6ris6 en ce que les donn^es^ g6om§triques 
comprennent des donnfees de vecteurs ( Cf - Cj> ^ ) entre deux 
noeuds donn6s. 

6. Proc6d4 selon I'une des revendications 1^5, 
25 caract6ris6 en ce que les donn^es g6om6triques 

comprennent des donn6es de vecteur normal (V) ^ la 
surface 4 nodSliser. 

7. Proc6d6 selon I'une des revendications 4^6, 
caract6ris6 en ce qu'au raoins certaines desdites donn6es 

3Q g^oB^triques sont assorties d'un coefficient ("or^^ . "CCT^p 
repr^sentatif du degr^ de certitude avec lequel ces 
donn6es sont connues et en ce que I'indice de violation 
des donn6es int6gre lesdits coefficients de certitude. 

8. Proc6d6 selon I'une des revendications 
35 pr^c^dentes, caract6ris6 en ce qu'il comprend en outre 

les Stapes consistant ^ mesurer des propri6t6s du corps h 



20 



FEUILLE DE REM PLACEMENT 



wo 91/04544 PCr/FR90/00676 

- 55 - 

trois dimensions dans la region d*au moins certains 
noeuds du maillage et ^ m^moriser les mesures obtenues 

(info) aux adressea sp6cifiques aux noeuds consid6r6s. 

9- Proc6d6 selon I'une des revendications 
^ pr6c6dentes , caract6ris6 en ce que, dans chaque 
coefficient de rugosity local, le rapport entre le poids 
(v« (k) ; a=k) du noeud concern^ et le poids (v^ (k) ; ae 
d'un satellite quelconque dudit noeud est 6gal au 
nombre (nw ) de satellites dudit noeud, affects d'un signe 
n6gatif. 

10. Proc6d6 selon la revendication 9, caract6ris6 
en ce que ladite combinaison pond6r6e utilise les memes 
poids • 

11. Proc6d6 selon I'une quelconque des 
revendications pr§c6dentes, caract6ris6 en ce que la 
somme des coefficients de rugosit6 locaux est une somme 
pond6r6e . 

12. Proc6d6 selon la revendication 11, rattach^e k 
la revendication 10, caract6ris6 en ce que ladite 

20 combinaison pond6r6e utilise 6galement les poids (ja(k)) 
associ6s aux coefficients de rugosit6 locaux. 

13. Proc6d6 selon I'une des revendications 
pr6c6de^tes., caract6ris6 en ce que la surface est une 
surface der separation entre deux milieux de propri^t^s ou 

25 de natures diff6rentes dans une formation g6ologique et 
en ce qu'on utilise la rapr6sentation de la surface pour 
optimiser 1 ' exploration de richesses souterraines de 
ladite formation. 

14. Proc6d6 selon I'une des revendications 1 A 12, 
caract6ris6 en ce que la surface est representative des 

^° variations d'un propri6t6 d'un corps k deux dimensions au 
voisinage d'un plan traversant ledit corps, ladite 
propri6t6 6tant representee par la dimension restante, et 
en ce qu'on utilise la representation de la surface pour 
optimiser 1 ' exploration ou 1 ' exploitation ' de r i chesses 
souterraines de ladite formation. 
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15. Proc6d6 selon I'une des revendications 1 ^ 14, 
caract6ris6 en ce que la representation de la surface est 
une representation graphique sur un support plan. 

16. Proc6d6 selon 1 ' une des revendications 1 4 14 
5 caract6ris6 en ce que la representation de la surface est 

une maquette en trois dimensions. 

17. Dispositif pour I'obtention d'une modeiisation 
d'une surface (S) representative notamment de 1' interface 
entre deux zones de natures ou de propriet6s differentes 
dans un corps ^ trois dimensions tel qu'une formation 
geologique ou un corps d'etre vivant, caracterise en ce 
qu'il comprend : 

- des moyens de mesure un ensemble de donn6es 
geometrigues relatives A la surface et associ6es 
respect ivement k des points donn6s de ladite surface; 

- des moyens pour effectuer un maillage de la 
surface de telle sorte que 1' ensemble desdits points 
donnes soit un sous-ensemble des noeuds du maillage; 

- une memoire contenant ^ des adresses donnees , 
20 specif iques k chaque noeud du maillage. les informations 

suivantes: 

* les coordonn6efi (Cff, C(^,cf|) du noeud consid6re, 

* le nombre (nb-sat) de noeuds satellites du 
noeud considere, 

25 * information (sat) permettant d'acc6der aux 

adresses sp6cifiques desdits noeuds satellites et par 
suite aux informations relatives k ces derniers, 

* le cas ech6ant des donnees geom6triques 
(const) associees audit noeud considere, 

2^ - des moyens de calcul pour ajuster de facon 

iterative les coordonnees de chacun des noeuds dont on ne 
connalt pas les coordonn6es pr6cises, en utilisant k 
chaque fois 1' addition d'une combinaison ponderee _des 
coordonnees actuelles des satellites et des satellites 
des satellites dudit noeud et d'une combinaison des 
donnees g6oin6tr iques associ6es audit noeud, de maniere k 
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amener k un minimum une somme (R*(Cp)) d'un indice de 
rugosity global (R(C^)), constitu6 par I'addition 
d'indices de rugosity locaux (RCCpjk)) respectivement 
associ6s k chacun des noeuds et chacun obtenu k partir 

5 d'une somme pond6r6e des coordonn6es actuelles dudit 
noeud et de ses satellites, et d'un indice global 
de violation desdites d6nn6es g6om6triques ; et 

- des moyens pour cr6er une repr6sentation de la 
surface k partir des coordonn6es ajust6es de chacun des 

10 noeuds • 
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